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I. INTRODUCTION 



The theory of quantum open systems has attracted 
growing attention because the framework of system plus 
environment has become a fundamental paradigm in 
quantum dissipative system, quantum optics and quan- 
tum information science [11, 0]- In reality, no quantum 
systems can be completely isolated from their environ- 
ment. The mutual interaction between the system and 
environment is responsible for many important physical 
processes such as dissipation, fluctuation, as well as deco- 
herence and disentanglement. Quantum open system is 
also an essential part of quantum information processing 
where the quantum noise and control are major issues 
in realizing quantum computing, quantum teleportation 
and quantum cryptography Q. 

While the theory of open systems in the Markov 
regimes can be extensively treated with the standard 
Lindblad master equations ll| or the corresponding 
Markov quantum trajectories |^-ul , the microscopic non- 
Markovian theory has not been well developed for the 
systems consisting of interacting spins, which are of ex- 
perimental interest Q. The difficulty arising from the 
strong coupling between the system and its environment 
is that the non-Markovian effects due to action and back 
reaction between the system and the environment must 
be taken into account. Undoubtedly, some memory re- 
lated properties such as quantum entanglement, inhib- 
ited spontaneous emission, quantum transport, quantum 
tunneling [9l-[l7| which are of importance will be lost 
when simply taking the Markov limit. Among all of 
the techniques attempt to investigate the non-Markovian 
evolution, the quantum state diffusion (QSD) equation 
initially proposed by Diosi, Gisin, and Strunz |18j pro- 



vides a powerful tool to deal with the quantum open sys- 
tems coupled to a boson bath [T9l - [2l] |. 

The recent developments of the non-Markovian dy- 
namics are concerned with a theoretical derivation of the 
QSD equation directly from an underlying microscopic 
model irrespective of environmental memory time and 
coupling strength 18|. The non-Markovian QSD equa- 
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tion is equivalent to the master equation formalism, but 
offers numerical advantages over the corresponding mas- 
ter equation if it exists. Typical to non-Markovian evo- 
lution is that the dynamics is explicitly dependent on 
the past history manifested in the integrals over the time 
in the dynamic equation [22h24j . In the case of non- 
Markovian QSD, under certain conditions^the nonlocal 
integral can be cast into a time-local form [isl. [2TI. [25l [26| . 
Then, the time-local QSD equation can be efficiently im- 
plemented as an analytical and numerical tools to de- 
scribe the non-Markovian effects. 

The primary purpose of this paper is to study the non- 
Markovian quantum dynamics of interacting qubits. We 
derive a closed exact time-local QSD equation for two 
interacting qubits coupled to a bosonic environment by 
replacing the functional derivative in the QSD equation 
with a time-dependent operator called O operator. We 
have established the explicit equations of motion for the 
coefficient functions contained in the O operator. The 
second purpose of this paper is to investigate the non- 
Markovian dynamics of qubit systems by using the de- 
rived exact QSD equation. Since the exact QSD equation 
is valid for an arbitrary structured boson bath and the 
strong coupling regime, we are thus able to address the 
important issue of entanglement generation via environ- 
ment memory in the full non-Markovian regime. It is 
shown that several interesting features may arise from 
this coupled qubit system when we examine the dynam- 
ics in a non-Markovian regime. With this coupled qubit 
model, we have studied the approximation based on the 
functional expansion of the O operator, and reveal some 
interesting features on non-Markovian perturbation the- 
ory. 
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The organization of this paper is as follows. In Sec. II, 
we introduce an interacting qubit model and develop the 
theory of the exact time-local quantum state diffusion 
equation. In Sec. Ill, based on the exact QSD equa- 
tion derived, we consider memory-assisted entanglement 
to illustrate the implications of non-Markovian dynam- 
ics on entanglement generation. As an example, we show 
that entanglement generation may be significantly mod- 
ified by the environment memory time. With the quan- 
tum trajectories generated by the exact QSD equation, 
we also analyze the steady state solution for the single 
trajectories recovering the steady state entanglement de- 
scribed by the density matrix. We conclude in Sec. IV. 

II. THE MODEL AND TIME-LOCAL 
QUANTUM STATE DIFFUSION EQUATION 

The model consists of a pair of coupled two-level atoms 
(or spins) interacting with a common bath. The total 
Hamiltonian can be written as (setting h — I), 

Htot = Hsys + Hcnv + Hint, (1) 

where 

Hsys = ^Acr^ + ^Bcrf + Jxyicricr^ + (^^(7+) + Jz<J^<^z , 
Hcnv = y^^ujjb]bj, 

3 

Hint = Y.{9,b]L + g*h,L^l (2) 

3 

where L — kao'^ + hb<^- is the system Lindblad oper- 
ator coupled to the environment. Note that ka,hb are 
constants describing different coupling strengths for two 
spins. 

This type of interaction between qubits has been stud- 
ied widely in the Heisenberg spin chain models [27l - [29| 
and proved to be a useful model in quantum information 
processing based on spin chain. Moreover, some cavity- 
QED systems may be described by a similar type of 
Hamiltonian in simulating the Heisenberg spin chain 'so'] . 
Notably, our model in the framework of quantum open 
system is a very general case, including many other phys- 
ically interesting models as special cases of this Hamilto- 
nian. For example, when Jz = 0, this model will reduce 
to the example of the Heisenberg XX type interaction. 
When Jxy = Jz = 0, two qubits will be treated as non- 
interacting spin. When lub = Jxy — Jz — k,b — 0, this 
model simply reduces to the one qubit dissipative model 

The stochastic Schrodinger equation for this interact- 
ing spin model may be formally written as fl§L [l9l [2l| , 

^ipt = -iHsysipt + Lzf-tjjt - J dsa{t, s)j^, (3) 

where a(t, s) — ISiP"^ is the bath correla- 

tion function and = —iJ2j9j^j^^'^^'^ is a complex 



Gaussian process satisfying M[zt] = M[ztZs] = 0, and 
M[ZjZs] = a{t,s). Here M[-] denotes the ensemble av- 
erage over the classical noise Zf. The quantum trajec- 
tory ^t{z*) recovers the density operator of the system 
in the ensemble average: pt — M[\ipt{z*)){ipt{z*)\] = 

The above exact equation contains a time-nonlocal 
term, it renders practical application impossible with- 
out evoking certain approximation. The existence of 
the time-local equation depends on if one can replace 
the functional derivative in the integral in Eq. ([3]) 
with a time-dependent (even noise-dependent) operator 

oit,s,z*y. 

— = 0{t,s,z )i;t{z ). (4) 
Using the "consistency condition" [Tsj 

5z* dt ~ dt Sz* ' ^ ' 

one can obtain a formal evolution equation for the oper- 
ator 0(t, s, z*): 

-0=[~iHsys + Lz;-L^O,0]-L^—0, (6) 

where 0{t,z*) = / dsa{t — s)0{t, s, z*). This equation 
Jo 

of motion for the O operator has to be solved with the 
initial condition, 

0(t,s = t,z*) ^ L. (7) 

The explicit O operator in Eg. Q has been determined 
for several interesting models [iMlUllHI. However, the 
construction of the O operator for the coupled two spins 
has not been found before. It is generally very difficult 
to find the exact O operator. As the major result of 
this paper, we have found the exact O operator for this 
coupled qubit model. To begin with, we note that the O 
operator may be expaned in the following way, 

Oit, s, z*) = hit, s)0, + /2(i, S)02 + hit, S)03 

+ hit,s)Oi+if ds'hit,s,s')z:,05, (8) 
Jo 

where the basis operators are given by 

Oi^afa^, 05^2a^a^, (9) 

and fj ij = 1,2,3,4,5) are some time-dependent coeffi- 
cients. Substituting Eq. ^ into Eq. we obtain the 
partial differential equations governing the coefficients of 
the O operator. 
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+2iujAfi + KAFifi - 


f hbF^/i - 


- iJxyh - KsFifa H 


- KsFifs 


+ 2a, /4 


+ KAFifi + 


K-sFsfi - 


ikbF^, 


|/.(M) = 


+ 2iuJBf2 + HAFif2 - 


f KBF2f2 - 


f 2i J,/3 + - 


\- kbFs/s 


iJxyfi 


- KAF2fi + 


KAF3.fi - 


ikaF^, 


|/3(M) = 


-iJxyfl - KA-F2/1 4 


- kaFs/i 4 


' 2lJj2 + KAFif2 + 




- 2iujBf3 


+ KAFifs + 


KBF2f3 - 


ikaF^, 




+2iJJi + KAFifi 4 


- KsF^fi - 


- iJxyf2 - 4 






+ KAFifi + 


KBFsfi ~ 


ikbF^, 


f5it,S,s') = 


+^.4^5/1 + 


- i^bFs/s 


- kaF^/a + "^ii^Afb 




4- KAFifi + KAFifb 


+ KBF2f5 


+ KBF3 



where Fj(t) — dsa(t, s)fj{t, s) (j = 1,2,3,4) and 

^5(^,5') — dsa{t, s)f5{t, s, s'), with the mitial condi- 
tions for Eq. PU)) : 

fl{t,S = i) = KA, 
f2{t,S = t) = KB, 
f3it,S = t) = 0, 

fi{t,s = t) = 0, 
/5(t,s-t,s') -0, 

fsit, S, s' ^t) = -i[KAf3{t,s) + KBfiit, s)]. (11) 

The time evolution for this model can be solved numeri- 
cally with the coefficients determined by Eq. (fTO| . 

With the explicit O operator the linear QSD equa- 
tion can be written compactly as a time-local equation: 



^ / 



(12) 



where 
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0{t, z*) = V Fj{t)Oj + i / ds'Fs 



{t,s')z:,0,. (13) 



Although the time evolution is, in principle, governed 
by the linear QSD equation ([3]), it is much more effi- 
cient to use the nonlinear version in numerical simula- 
tions since it preserves the norm of the state vector. The 
nonlinear QSD equation can be written as [Tsj . 



d_ 

dt 



- [(Lt - (Lt) JO - ((Lt - {L^)^)0)]^u (14) 



where ipt = is the normalized state vector , (L)^ = 

ipt L Tpt^ denotes the mean value of the operator and 

Zf — 2j -I- Jg dsa*{t, s) {L^)g is the shifted noise. 

Let us emphasize that our results are valid for an ar- 
bitrary correlation functions. For simplicity, in this pa- 
per, we model the noise as an Ornstein-Uhlenbeck process 



for which the correlation function is a{t,s) = ^e"''''*"*'. 
This correlation function is very useful in demonstrating 
the non-Markovian effects as well as the Markov case by 
varying the parameter 7, which controls the memory ef- 
fect of the reservoir. When 7 — )■ 00, a{t, s) —>■ 6{t — s) it 
corresponds to the Markov case. Oppositely, when 7 is 
small, the system is in a non-Markovian regime. 



III. NUMERICAL RESULTS AND 
DISCUSSIONS 

Having developed the time-local non-Markovian QSD 
equation for the two interacting qubits, this section is 
concerned with the non-Markovian dynamics of entan- 
glement involving entanglement generation and evolution 
[3ll - [39l l. Throughout the paper, the degree of entangle- 
ment is measured by "concurrence" [ioj . 



A. Memory-assisted entanglement generation 



With our exact non-Markovian QSD equation (|T4|) . we 
are able to consider the important non-Markovian fea- 
ture of entanglement generation. For this purpose, we 
consider the separable state tpo = |11) to show how the 
environmental memory affects the entanglement genera- 
tion. It has been known that the entanglement may be 
generated for a qubit system through coupling qubits to 
a common bath. However, it is not clear how the environ- 
mental correlation times affects entanglement generation. 
In Fig. [U we plot the entanglement evolution for differ- 
ent parameters of 7, where I/7 gives the memory time of 
the environment. With this initial state, we can see that 
when 7 > 2 (close to the Markov limit for this model), 
there will be no visible entanglement generated for all 
times. In fact, it is easy to show that in the Markov 
limit no entanglement will be generated (i.e. 7 — >■ 00) 
for this special initial state. In addition, it is seen that 
when the parameter is too small, the entanglement gen- 
eration becomes less significant shown in our numerical 
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FIG. 1: Entanglement generation for different memory times 
of the noise. A particular initial state is chosen as ipo = • 
The curves for different 7 are marked in the graph. The other 
parameters are lja = ujb ~ 0.5tJ, ka ~ kb ~ 1, Jxy = 0.5. 
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FIG. 2: Entanglement generation vs 7 and t. A particular 
initial state is chosen as t/jq = . The other parameters are 

UJA = IJJB — O.buj, Kl = HB — 1, Jxy = Jz = 0. 



simulations. This shows that the non-Markovian prop- 
erty affects generating entanglement in a subtle way. 

We can show that entanglement generation for this 
special initial state is due to the environmental memory 
rather than the spin-spin coupling. Interestingly enough, 
we find that the largest amount of entanglement may be 
generated if the correlation time of environment is nei- 
ther too long (7 is small) nor too short (7 is large). In 
Fig. [21 we plot the case without coupling between the 
two spins. The results confirm that entanglement for the 
two non-interacting qubits is purely generated by envi- 
ronment memory. Our results here have clearly shown 
that non-Markovian noise can be very useful rather than 
harmful to the entanglement if the memory times are 
properly chosen. 



FIG. 3: Entanglement generation. The red (solid), green 
(dashed) and blue (short-dashed) lines represents 7 = 0.5, 
7 = 1, 7 = 2, respectively. The initial state is tpo = |10) . The 
other parameters are uua = <^b ~ 0.5aj, Jxy = 0.5, Jz ~ 0.1, 
KA = HB — 1. 



B. Steady state entanglement 

Quantum trajectories generated by the non-Markovian 
QSD equation ([3]) have many interesting properties. For 
instance, we can use those quantum trajectories to study 
the steady state entanglement. In Fig. [3J we plot the 
entanglement dynamics from a separable initial state 
4>Q — 1 10). The numerical results show that the initial 
pure state evolves into a steady entangled state. In fact, 
with the help of non-Markovian trajectories, we can show 
that all the trajectories will localize to the eigenstates of 
the Lindblad operator L and the residual entanglement 
of the steady state is not zero. Now, let us examine 
the large t behavior of the single trajectories. First, we 
choose the following four states 

l'/'i) = |ll), 

1^-2) = ^(|10) + 101)), 

|V'3) = ^(|io)-|oi)), 

IV4) = |00) . (15) 

which form a complete basis of the Hilbert space. Hence, 
an arbitrary state can always be written as 

IV-) = Ci \^Pi) + C2 \lp2) + C3 IV-S) + C4 IV'4) , (16) 

where jcip + |c2p + jcsp + |c4|2 = 1. 

Substituting Eq. ^ into the QSD Eq. dH]), we get 
the following stochastic differential equations for the co- 
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efHcients: 

ot 
d 

— C2(i) = i{Jxy - Jz)c2 + \/2z(*ci - 2Fi{t)c2 

z*)ci, 

d 

— C3(i) = -l{Jxv + Jz)c3, 

d 

— C4(i) = i{Jz - 2uj)c4 - V2Z( C2. 



(17) 
(18) 

(19) 
(20) 



where F^{t,z*) — i j^ds' F^(t,s')z*,. Here, we assume 
that the two qubits are identical and symmetric {uia = 
^BtI^a = i^b) for simphcity. The more general case can 
be dealt with in a similar way. For the initial state tp^ = 
1 10), the initial conditions for this set of equations are 



For the above initial conditions, ci{t) = at any times, 
and the steady state solution for C2 {t) is simply given by 



C2(oo) = 0. 



(25) 



Thus, the large t state for a single trajectory will not 
contain |-!/)i) and 1-02) components, i.e. 



I^final) = C3 IV's) + C4 |?/'4) 



(26) 



Now, it is clear that all trajectories will always evolve into 
a sub-space of the Hilbert space spanned by j-j/'s) and 
which are the eigenstates of the Lindblad operator L = 
CT^ + ct:^. Moreover, we can easily derive the equations 
of the motion for the ensemble means of the cross terms. 
For example, M[c3C4] satisfies the following equation: 



di 



M[c3C4] = {iJxy + '2iJz - 2iw)Af [C3C4] - V2M[z*clc2] 



Since 02(00) = in the long-time limit, we have 

M[c3C4] -> 0, A/[c3C4] ^ 0. (27) 

where t is large. Therefore, the density matrix recovered 
by pt = M[\^pt{z*)){'ipt{z*)\] will approach to the steady 
state density matrix in long-time limit: 







M[4c3]\t^ 










Af[c|C4] 



(28) 



in the basis {iV'i) , IV'2) , IV's) j IV'4)}- The steady state 
density matrix in the basis {|11) , |10) , |01) , |00)} is then 




FIG. 4: Time evolution of the fluctuation (AL)^ and the 
probability |c3p -|- |c4p for quantum trajectories. Each curve 
represents a single quantum trajectory. The parameters are 

UJA — iJ^B = O.SOJ, KA = K.B = 1, "t — 1, Jxy — 0.5. 



ci(0) 


-0, 


(21) 


C2(0) 


1 


(22) given by 


C3(0) 


1 

"71' 


(23) 


C4(0) 


= 0. 


(24) 





P22 P23 

Ph P33 

P44 



(29) 



where p22 = P33 = ^Af^sCs] |t->oo, P23 = 
~^M[clc3]\t^oo and ^44 = M [c|c4]|t^oo■ 
The analytical results shown above can be verified by 
our numerical simulation as illustrated in the figure 4. 

In Fig. |4] (a), we plot the fluctuation of the system 
Lindblad operator L defined as (AL)^ = (L-^) — (i) . It is 
shown that, as time evolves, the fluctuation of L converge 
to zero in the final steady state for all trajectories. In Fig. 
|4](b), we plot the sum of the probabilities of ji/'a) and \ip4) 
for a few trajectories. The sum of the probability |c3p -I- 
|c4p will always evolve to 1 for any single trajectories 
which means ci and C2 are always zero in the final steady 
state. The numerical results in Fig. |4] simply show that 
the final state of single trajectories take the form of Eq. 
(|26p as predicted by analytical analysis. In addition, our 
numerical results show that the final density matrix is 
approximately given by 



Pfinal 




0.25 
-0.25 






-0.25 

0.25 

0.5 



(30) 



which is also consistent with our prediction in Eq. 

It should be noted that such a steady state has be 
shown to be the stationary solution of the Markov mas- 
ter equation (41^^^]. In our QSD method, it is showed 
that this state is also the stationary solution of the non- 
Markovian quantum trajectories. Therefore, the ensem- 
ble mean of the quantum trajectories recovers the non- 
Markovian steady state of the density matrix. It is im- 
portant to note that our non-Markovian theory does not 
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FIG. 5: Time evolution of purity for different initial states. 
The initial states for red (solid), green (dashed), blue (short- 
dashed) and black (dash-dotted) lines are ^pa = --1=(|11) + 

|00)), = ;^(|01) + |10)), i^o = |10), and i/>o = 
respectively. The other parameters are uja ~ (j^b ~ 0.5cj, 
KA = KB = 1, 7 = li Jxy = 0.7, Jz = 0.3. 

depend on the existence of the non-Markovian master 
equation. 

Finally, we present the time evolution of the purity, de- 
fined as Purity ~ tr{p'^), vifhich determines whether the 
quantum state is pure or mixed. Specifically, we find the 
purity of the initial state ij^Q = |10) i which evolves into 
a steady state with some residue entanglement preserved 
in the final steady state. In Fig. [5l it is validated that 
the final steady state is a mixed entangled state since 
the purity is smaller than one at the steady state. The 
purity is approximately 0.5 in long time limit, which is 
consistent with Eq. ((30|). For all the other three initial 
states, we will see two qubits will firstly evolve into a 
mixed (perhaps entangled) state by the coupling interac- 
tion, then decay into the final ground state |00) (which 
is a pure state) due to the dissipative environment. 

C. Exact versus approximate QSD equations 

For the coupled two-qubit model, the exact O operator 
contains five terms, as shown in Eq. ([8]). It is interesting 
to investigate the perturbative QSD equations based on 
the exact expression of the O operator. For example, we 
may drop the first order noise-dependent term O5 and 
only keep the first four terms, i.e. the O operator in Eq. 
([5]) is approximately written as 

0(i,S,Z*)«/l(t,s)Ol+/2(t,s)02 

+ h{t,s)Os + h(t,s)Oi. (31) 

Clearly, the first order term O5 becomes less important 
when the system is close to the Markov regimes. Also, nu- 
merical calculations will be greatly simplified if the first 
order noise is neglected. In Fig. [6l we have compared the 
exact time evolution of entanglement and the approxi- 
mate case without O5. The numerical result shows that 




tot mt 



FIG. 6: Comparison of the dynamics generated by the ap- 
proximate O operator and the exact O operator. The red 
(solid) curve represents the exact QSD equation, and the 
green (dashed) curve is for the approximate QSD. (a) (b) (c) 
and (d) represent four different initial states -^(llO) 4- |01)), 
-^(|11) + |00)), 1 11) and |10), respectively. The parameters 
are chosen as 7 = 0.3, oja = i^s = 0.5tJ, = hb ~ 1, 

Jxy ■— 0.5. 



the accuracy of the zeroth order approximation is sensi- 
tively dependent on the initial state of the qubits. For 
the anti-correlated initial states. Fig. |6](a) and (d) shows 
that the approximation is extremely good compared with 
the exact result. On the other hand, for the correlated 
initial states in case (b) and (c), the noise dependent 
term tends to smooth out the entanglement curve. In 
fact, this can be explained by analyzing the QSD Eq. 
(|17ll20p . If the initial states do not contain ipi compo- 
nent, i.e. ci(0) = ('0o|'0i) — Oi then ci{t) is always zero 
during the time evolution. The fifth term of the O op- 
erator O5 = 2(7^(7^ will not affect the dynamics of the 

system, since the coefficient of O5 , i.e. i^5(t, z*), is al- 
ways associated with ci(t). This observation explains our 
results plotted in Fig. 6. As shown in Fig. 6, the noiseless 
approximation (setting O5 = 0) gives extremely good re- 
sults for the two initial states which do not contain |11) 
component. Our simple analysis here has already shown 
the complexity level of the approximate QSD equations 
arising from the noise expansion of the O operator. More 
in-depth study on the approximate QSD equations will 
be presented elsewhere. 
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IV. CONCLUSION 

In conclusion, the central aim of this paper is to de- 
velop a time-local theory of non-Markovian quantum 
state diffusion for a model consisting of two interacting 
qubits - a model that is of interest for both quantum 
information processing and quantum optics. As an ap- 
plication, we used the non-Markovian quantum trajec- 
tory to analyze the entanglement evolution of the qubit 
system for several parameters ranging from strongly 
non-Markovian to Markov regimes. Our results also 
showed explicitly how the environmental memory gen- 
erates the entanglement of the qubit system. With the 
non-Markovian QSD equation, we further studied the en- 
tanglement evolution and the steady entangled state for 



the QSD equation. Our time-local QSD approach may 
be extended to multiple qubit systems, this will be the 
topic for a future publication. 
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